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Abstract 



The strange quark mass is determined from a new QCD Finite Energy Sum 
Rule (FESR) optimized to reduce considerably the systematic uncertainties 
arising from the hadronic resonance sector. As a result, the main uncer- 
tainty in this determination is due to the value of Aqcd- The correlator of 
axial-vector divergences is used in perturbative QCD to five-loop order, in- 
cluding quark and gluon condensate contributions, in the framework of both 
Fixed Order (FOPT), and Contour Improved Perturbation Theory (CIPT). 
The latter exhibits very good convergence, leading to a remarkably stable 
result in the very wide range Sq = 1.0 — 4.0 GeV^, where Sq is the ra- 
dius of the integration contour in the complex energy (squared) plane. The 
value of the strange quark mass in this framework at a scale of 2 GeV is 
m,(2 GeV) = 95 ± 5 (111 ± 6) MeV for Aqcd = 420 (330) MeV, respectively 
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1 Introduction 



The strange quark mass is a very important QCD parameter measuring the 
strength of chiral SU{3) ® SU{3), and flavour SU{3) symmetry breaking. It 
also has a strong impact on a variety of QCD applications in weak hadronic 
physics. For this reason, many attempts have been made in the past to de- 
termine rris in various frameworks, e.g. QCD sum rules [I]- [2], and Lattice 
QCD [3j. The most recent QCD sum rule determinations from the pseu- 
doscalar channel have made use of state of the art results in perturbative 
QCD (PQCD) to five-loop order [1]. In spite of this, the real uncertainty in 
the value of remains high due to the hadronic resonance sector. In fact, 
beyond the kaon pole, the pseudoscalar hadronic spectral function is not 
known from direct experimental mesurements. Two radial excitations of the 
kaon have been observed [5], and a certain amount of theoretical input has 
gone into attempts to build a reasonable spectral function incorporating these 
resonances. However, inelasticity and non-resonant background are realisti- 
cally impossible to model. This constitutes a form of systematic uncertainty 
seriously limiting the precision of these determinations. In summary, current 
information on the QCD side of the sum rules is not matched in quality by 
the pseudoscalar hadronic sector. An attempt to rectify this situation has 
been made recently [6j in the form of a new QCD Finite Energy Sum Rule 
(FESR) involving as integration kernel a second degree polynomial which is 
required to vanish at the peaks of the two pseudoscalar resonances. As a 
result of this, the kaon pole and the QCD contributions dominate the FESR; 
the importance of the hadronic resonance sector being reduced by up to an 
order of magnitude. This FESR was used in [B] to determine the scalar and 
pseudoscalar correlators at zero momentum, and the strange quark conden- 
sate. An upper bound on the strange quark mass was also obtained there, 
e.g. for the running mass at a scale of 2 GeV this bound is 



In this paper we use this FESR in the pseudoscalar channel to determine the 
value of strange quark mass to five-loop order in PQCD, and including the 




121 MeV (Aqcd = 330 MeV) 
105 MeV {Aqcd = 420 MeV) . 
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leading vacuum condensates. We use the framework of Fixed Order Pertur- 
bation Theory (FOPT), as well as Contour Improved Perturbation Theory 
(CIPT). We find that as a result of the integration kernel in the FESR, the 
hadronic resonance contribution is considerably reduced relative to the kaon 
pole. The latter is of the same order as the PQCD contribution, and both are 
up to a factor five bigger than the resonance term, this being comparable to 
the gluon condensate term. Numerically, the resonances add roughly 10% to 
the quark mass, relative to the value obtained from the kaon pole and PQCD, 
while the gluon condensate subtracts a similar amount, and the light-quark 
condensate reduces it by another 1-2%. Higher dimensional condensates and 
higher order quark-mass terms contribute negligible amounts. Results from 
FOPT for the running strange quark mass at a fixed scale are reasonably 
stable in a wide range of values of sq, the upper limit of integration in the 
FESR (sq ~ 2.5 — 4.5 GeV^). However, in CIPT the stability is remarkable 
in an even wider range, e.g. ms{2 GeV) changes by less than 1% in the 
range Sq ~ 1.0 — 4.0 GeV^. There is a very strong correlation between 
and the value of the QCD scale A, which produces most of the uncertainty 
in the result for the strange quark mass (roughly 16 %). However, unlike 
the situation in the hadronic resonance sector, the uncertainty in A can be 
reduced, in principle, by improving the accuracy of the theoretical input, as 
well as of the data used in its determination. 



2 Fixed Order Perturbation Theory 

We first introduce the correlator of axial-vector divergences 

Mq') J d'x e^'^^ < \T{d^A,ix) , 9^4(0))| > , (2) 

where d^^Afj_{x) = {vHs + m„) : 's{x)i^r^u{x) : is the divergence of the axial- 
vector current. To simplify the notation we shall use in the sequel rris+mu = 
m. Finite Energy Sum Rules (FESR) involving this correlator follow from 
Cauchy's theorem in the complex energy-squared, s - plane (see Fig. 1), i.e. 
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and the contour integral is performed over a large circle where the exact 
-05(5) can be safely replaced by its QCD counterpart ip^^^ls). We introduce 
now an integration kernel in the form of a second degree polynomial 



A5(s) = 1 - ao s - ai , (4) 

where ao, and ai are free parameters to be fixed by the requirement that 
A5(M^) = A5 (M|) = 0, with Mi^2 the masses of the two resonances in the 
strange pseudoscalar channel, K(1460) and K(1830) [5]. This gives 



ao = 0.768 GeV-2 = -0.140 GeV"^ . (5) 
Writing explicitly the kaon pole, Cauchy's theorem, Eq. 3, becomes 



dsij^^'''{s)A,{s) = 2flMt,A,{Ml) 



2m 



C{\so\) 
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+ I (is - /m ^5(5)1 ijij5 As (s) , (6) 



where Sth = [M^ + M^^y is the resonance threshold, and fx = (1.22 ± 
0.01)/^ = 113 ± 1 MeV. The introduction of the integration kernel A5(s) is 
expected to reduce the importance of the hadronic resonances in the deter- 
mination of rris. A posteriori, this is in fact the case: in the stability region 
for ms(so) the resonance contribution is roughly a factor five smaller than 
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Figure 1: Integration contour in the complex s-plane. 



the pole term, as well as the PQCD piece. 

In the framework of FOPT as and rus are taken as constants, and only 
terms involving powers of log{—s/fi^) contribute to the contour integral. The 
renormalization group summation of leading logs is only carried out after the 
contour integration by setting //^ = —sq. The PQCD result for ipbi^) up to 
four- loop order has been known for quite some time [7], while the five-loop 
expression for its second derivative has been obtained recently [1]. Integrating 
the latter twice provides ip^i^s) up to a non-contributing polynomial. We 
define 




(7) 



and obtain in PQCD the following results 
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where {^G"^) ~ 0.06 GeV^ and (gg) ~ (-250 MeV)^ [9]. We find that 
terms of 0{m^) and higher, as well as the strange quark condensate, give 
negligible contributions in the region of stability (which turns out to be 
So — 2.5 — 4.5 GeV^). This is also the case for the condensates of dimension- 
six and higher. 



Turning to the hadronic sector, the spectral function in the pseudoscalar 
channel involves in addition to the kaon pole, at least two radial excitations, 
the K(1460) and K(1830) both with widths of about 250 MeV [5]. We follow 
the procedure outlined in [10], where the resonance part of the spectral func- 
tion is written as a linear combination of two Breit-Wigner forms normalized 
at threshold according to chiral perturbation theory. The latter incorporates 
the resonant sub-channel i^*(892) — vr which is important due to the narrow 
width of the i^*(892). Other embellishments are certainly possible, but then 
the presence of the integration kernel A5(s) in Eq.(6) makes these attempts 
unnecessary. In fact, the resonance contribution to Eq.(6) in a wide range of 
values of sq is up to a factor five smaller than the PQCD term, and a similar 
factor smaller than the kaon pole contribution. 



Adding up Eqs.(8)-(14) gives the left hand side of the FESR, Eq.(6), which 
leads to the results for ms(2 GeV) shown in Fig. 2. Curve (a) corresponds to 
A = 330 MeV, and curve (b) to A = 420 MeV. In obtaining the results shown 
in Fig. 2, the light quark mass has already been subtracted; its value can be 
safely estimated from the chiral perturbation theory ratio [TT]: ^ = 24.4 ± 
1.5, where = (m^ -|- mrf)/2. For the running coupling needed in Eqs.(9)- 
(14), enough precision is obtained using the four-loop result. We have also 
achieved enough precision with the four-loop result to convert the running 
quark mass m7(so) into ms(2 GeV). Numerically, the latter is ms(2 GeV) 
= 114±8 MeV (102±6 MeV), for A = 330 MeV (A = 420 MeV), respectively. 
The breakdown of the various factors contributing to the quark mass is as 
follows. To the basic result from PQCD and the kaon pole, the hadronic 
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Figure 2: Running strange quark mass in FOPT at a scale of 2 GeV as a 
function of Sq. Curve (a) is for A = 330 MeV, and curve (b) for A = 420 
MeV. 



resonances add roughly 10 % to the mass, the gluon condensate reduces the 
mass by a similar amount, and the light quark condensate by an additional 
1-2 %. The errors quoted above are only due to variations of ms{2 GeV) 
inside the stability region Sq = 2.5 — 4.5 GcV^. A more realistic error can be 
established by assuming that the resonance contribution to the quark mass 
has been understimated or overestimated by a certain amount. Choosing 
this amount conservatively as 30 % produces an additional uncertainty in 
rus (2 GeV) of ±4 MeV. In addition, assuming that the unknown six-loop 
contribution is comparable to the five-loop one would reduce (2 GeV) 
by roughly 2 MeV. Combining all uncertainties and using the breakdown 
described above, the final results for the quark mass are 
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m,(2 GeV) 



114 ± 14 MeV (Aqcd = 330 MeV) 
102 ± 11 MeV (Aqcd = 420 MeV) . 



(15) 



3 Contour Improved Perturbation Theory 

Contour Improved Perturbation Theory [12j has been shown to provide better 
convergence than FOPT in the QCD analysis of the vector and axial-vector 
correlators in tau-lepton decay. We find this to be also the case for our 
analysis of the contour integral in Eq.(3). Unlike the case of FOPT, where 
as(so) is frozen in Cauchy's contour integral and the Renormalization Group 
(RG) is implemented after integration, in CIPT as is running and the RG is 
used before integrating. This is done through a single-step numerical contour 
integration and using as input the strong coupling obtained by solving nu- 
merically the Renormalization Group Equation for as{~s) . This technique 
achieves a partial resummation of the higher order logarithmic integrals, and 
improves the convergence of the PQCD series. CIPT has been used success- 
fully in QCD analyses of tau-lepton hadronic decays [9], [12]. In the case 
of the pseudoscalar correlator involving the running quark mass as an over- 
all multiplicative factor, implementation of CIPT requires that not only the 
running coupling but also the running quark mass be integrated around the 
Cauchy contour. The running quark mass can be computed at each step 
by solving numerically the corresponding RGE. To establish notation and 
conventions, we write the RGE for the coupling as 



where = Os/vr, and for three quark fiavours /?o = 9/4, /3i = 4, = 
3863/384, (3^ = (421797/54 + 3560C(3))/256. In the complex s-plane s = 
So with the angle x defined in the interval x G (— 7r,7r). The RGE then 
becomes 




(16) 



N=0 



9 



das{x) 



N+2 



(17) 



dx 



N=0 



This RGE can be solved numerically using e.g. a modified Euler method, 
providing as input as{x = 0) = as(— Sq)- Next, the RGE for the quark mass 
is given by 



where for three quark flavours 70 = 1, 71 = 182/48, 72 = [8885/9 — 
160 C(3)]/64, 73 = [2977517/162-148720 C(3) /27+2160 C(4)-8000 C(5) /3]/256. 
With the aid of Eqs. (16)-(17) the above equation can be converted into a 
differential equation for m{x) and integrated, with the result 



where the integration constant m(0) is identified as the overall multiplica- 
tive quark mass in the expression for the pseudoscalar correlator, i.e. m = 
[ms(so) +mu{so)]. Cauchy's theorem, and the resulting FESR will be written 
for the second derivative of ip5{s), in which case it is straightforward to show 
the following identity 




(18) 



M=0 




(19) 
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where 
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and g{s) is an arbitrary analytic function which we choose as g{s) = A^ls), 
with A^{s) given in Eq.(4). In this case instead of Eq.(6) the FESR becomes 
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with g2 = —q^, and Renormahzation Group improvement has been used 
to dispose of the logarithmic terms. The rather long four- and five-loop 
expressions can be found in |1] and [7j. The left hand side of Eq.(22) in 
PQCD can be written as 
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where Ko = Cqi, K, = Cn + 2Cu, = C^i + 2C22, = C31 + 2C32, 
= C41 + 2C42, with Cik defined after Eq.(12), and 
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and bi = —{sq — aosl/2 — aiSg/S), 62 = 1/2, 63 = — ao/6, and 64 = — ai/12. 
The contribution of the gluon condensate to the left hand side of Eq.(22) is 
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(28) 



where the scale /xq — 1 GeV^ appears in connection with the removal of 
logarithmic quark mass singularities (see [2]). The light-quark condensate 
contribution is given by 
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Figure 3: Running strange quark mass in CIPT at a scale of 2 GeV as a 
function of sq- Curve (a) is for A = 330 MeV, and curve (b) for A = 420 
MeV. 



Substituting Eqs. (26), (28) and (29) in the FESR, Eq.(22), completes the 
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expression giving the running quark mass m(so)- The breakdown of the 
various contributions is similar to the case of FOPT. The Kaon pole and 
PQCD are the leading terms in the FESR, the hadronic resonances add 8-9 
% to the mass, the gluon condensate reduces it by roughly the same amount, 
and the light-quark condensate induces a further reduction of 1 -2 MeV. 
As expected, the convergence of the PQCD series is much better than the 
FOPT counterpart. The results are shown in Fig. 3 for the running strange 
quark mass at 2 GeV. Curves (a) and (b) correspond to A = 330 (420) MeV, 
respectively. The stability of the results, and the wide range of the stability 
region are quite remarkable. In fact, (2 GeV) changes by not more than 
1-2 MeV in the interval sq — 1 — A GeV^. This stabihty is present at each 
loop level in PQCD, i.e. the quark mass as a function of Sq is essentially flat 
if computed at one loop level, one- plus two- loop level, etc.. 

To arrive at a reasonable error, we proceed as before in FOPT, and assume 
that the hadronic resonance contributions could be underestimated or over- 
estimated by up to 30 %. In this case, the uncertainty in rUg (2 GeV) for A = 
330 (420) MeV, would be ±3 (2) MeV. In addition, if one were to assume 
the extreme scenario in which the six-loop contribution would be compara- 
ble to the five-loop one, then a further reduction in the mass by roughly 3 
MeV would result for both values of A. Combining these two uncertainties 
linearly, the final results are 



A breakdown of the contribution of each loop-order in PQCD to the quark 
mass rus (2 GeV) is as follows. For A = 330 MeV, the pole plus resonances 
plus condensates and only the one-loop contribution would give (2 GeV) 
= 152 MeV, including the two-loop term reduces this to rris (2 GeV) = 129 
MeV, with the three-loop it gives (2 GeV) = 120 MeV, with up to four 
loops this becomes (2 GeV) = 115 MeV, and finally with all five loops 
the result is (2 GeV) = 111 MeV. Similar relative differences are found 
for A = 420 MeV. As seen from Eq. (30) the main uncertainty is due to 
the value of A. Finally, we have considered the impact on the above results 




111 ±6 MeV (Aqcd = 330 MeV) 
95 ± 5 MeV {Aqcd = 420 MeV) . 



(30) 
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from uncertainties in the vacuum condensates as well as inclusion of higher 
order quark mass corrections. Allowing the extreme case of an uncertainty 
of a factor two in the gluon condensate produces changes in ms{2 GeV) 
well within the errors given in Eq.(30). The contribution of the strange 
quark condensate is negligible, and that of the light (up- or down-) quark 
condensate, which is known accurately, is at the level of 1%. Higher order 
quark mass corrections, as well as vacuum condensates of dimension d = 6, 
have no impact on the results and can be safely ignored in this method. 
The results above satisfy comfortably the upper bounds given in Eq.(l). In 
this particular application, CIPT has proven to be far better than FOPT. 
Although the results from both methods agree within errors, CIPT leads to 
a remarkable stability of ms{2 GeV), in a remarkable wide range of values of 
So- For this reason, we would not advocate combining the results from both 
methods. 



4 Conclusions 



The main advantage of using pseudoscalar correlators to determine the quark 
masses is that they enter the QCD expressions as overall multiplicative fac- 
tors, rather than as corrections to a leading term. In addition, they in- 
volve the pseudoscalar pole with parameters well known from experiment. 
Unfortunately, there is no direct experimental information on the hadronic 
resonance spectral function, except for the masses and widths of the first 
few resonances. This information is not enough to reconstruct reliably these 
spectral functions, as inelasticity and non-resonant background effects are 
realistically impossible to guess. For this reason, quark mass determinations 
from pseudoscalar correlators are affected by endemic systematic uncertain- 
ties not subject to improvement. In this paper we have used a new QCD 
FESR [6j for the strange pseudoscalar correlator, involving as integration 
kernel a second degree polynomial which is required to vanish at the peaks of 
the first two radial excitations of the kaon. As a result of this, the relative im- 
portance of the hadronic resonance sector in the determination of the strange 
quark mass is considerably reduced. In fact, this contribution turns out to 
be up to a factor five smaller than the leading contributions from PQCD 
and from the well known kaon pole. We have used the techniques of FOPT 
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as well as CIPT to compute the integrals in the complex energy (squared) 
plane. The latter method is far superior to the former, and gives a running 
mass at a fixed scale which is remarkably stable in a very wide range of sq, 
the radius of the integration contour in the complex plane. For instance, for 
the strange quark mass at 2 GeV we find (2 GeV) = 111 (95) MeV, for A 
= 330 (420) MeV, respectively, in the wide range Sq = 1 — 4 GeV^. To arrive 
at a reasonable, but still conservative estimate of the uncertainties we have 
assumed that (a) the resonance parametrization might be an underestimate 
or an overestimate of the hadronic spectral function of up to 30 %, and (b) 
the unknown six-loop PQCD term could be comparable to the five-loop term. 
Each of these assumptions induces an uncertainty in the quark mass at the 
level of 3 %. Adding them linearly gives the results in Eq. (30). The main 
uncertainty is then due to A. However, unlike the hadronic resonance spec- 
tral function, this error is subject to improvement. The results from FOPT, 
Eq. (15), are expectedly in agreement with those from CIPT. However, in 
FOPT an additional sizable uncertainty arises from the variation of rris (2 
GeV) in the (narrower) stability range. Results from both methods satisfy 
the upper bound Eq. (1). Comparison of our results with previous determi- 
nations is made somewhat difficult due to various reasons. Some of the very 
old determinations were afflicted by logarithmic quark mass singularities in 
the correlators. This issue was only clarified in p]. In addition, the values of 
A used in the past were much lower than at present. Given the strong cor- 
relation between A and the quark mass, this becomes a serious issue. Next, 
knowledge of the PQCD contribution has improved considerably over the 
years, from two-loop level to the current five-loop level. Last, but not least, 
the systematic uncertainties due to a lack of direct experimental information 
on the hadronic resonance spectral function may have been underestimated 
in the past. In any case, comparing the results in Eq. (30) with the most 
recent determinations [1], [3j, [4j, shows very good overall agreement. 

In closing we wish to mention that in some applications of FESR, e.g. in 
tau-decay, perturbative QCD does not appear to hold close to the real axis. 
It is not entirely clear whether there is a problem with the data, or with 
PQCD itself. In any case this has led to the proposal of weighted FESR with 
weight functions vanishing at s = sq fL3j. These so called pinched FESR 
improve considerably the saturation of the Weinberg sum rules, and resolve 
some inconsistencies in the determination of vacuum condensates in the vec- 
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tor and axial-vector channels. 



With this background it is reasonable to investigate the impact of such ker- 
nels in the determination of the strange quark mass discussed here. First 
of all, the rate at which a hadronic spectral function approaches its PQCD 
limit is channel dependent. Presumably, potential duality violations share 
this feature. In the framework of the method discussed here, it turns out 
that the ratio of the QCD and the hadronic contributions, which equals the 
square of the running quark mass, leads to a value for ms(2 GeV) which is 
remarkably stable as a function of the upper limit of integration sq, in an 
also remarkable wide range Sq — 1 — 4 GcV^. Such quality is hardly found 
in typical QCD sum rule applications. It is then reasonable to conclude that 
in this particular channel, and using our integration kernel, duality appears 
to be well satisfied. It should be stressed that the motivation for introducing 
the kernel Eq.(4)in the pseudoscalar channel is rather different from that for 
the pinched kernel in tau-decays. In fact, the hadronic spectral functions in 
the latter case are known from experiment, while this is not the case for the 
pseudoscalar channel. 

In any case, and to continue looking at this issue, we may consider the di- 
rect product of the kernel Eq.(4) and a pinched one. The first undesirable 
result of such a procedure is that the kaon pole contribution now becomes 
a function of sq, and is numerically reduced. This does not happen to the 
Weinberg sum rules, as they are valid in the chiral limit, in which case the 
pseudoscalar meson pole contribution to the spectral function involves the 
delta function S{s). As a result of this, the pinched kernel does not affect this 
contribution. This behaviour of the kaon pole contribution is contrary to the 
spirit of the method used here. In fact, since there is accurate experimental 
information on this pole, one wishes to enhance its contribution rather than 
reduce it. In any case, using this additional kernel we have studied the con- 
vergence of the PQCD series, and the relative contributions of the various 
terms, e.g. vacuum condensates, higher order quark-mass corrections, pseu- 
doscalar meson pole and resonances. The result is that the addition of the 
pinched kernel has only a negative impact on the results. The convergence 
of the PQCD contributions is not as good, and the stability region is consid- 
erably reduced. Numerically, though, the change in the final value of ms(2 



17 



GeV) is well within the error given in Eq.(30), but with a narrower stability 
region. We must then conclude that there is no advantage in introducing an 
additional pinched kernel in this channel. 



Note added in proof: After this work was completed, a new and more accurate 
determination of ^^(Mr) by the ALEPH collaboration has been released [H], 
which implies a narrower range for A, i.e. A = 365 — 397 MeV. Using 
this in our determination we obtain ms{2GeV) = 99 ± 5(105 ± 6) MeV for 
A = 397(365) MeV, respectively. Combining both values gives ms{2GeV) = 
102 ± 8 MeV. 
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